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1 Introduction 

The aim of this paper is to present some classical and more recent results con- 
cerning the uniqueness of limit cycles for the Lienard equation 

x + fix)x + g{x)^0. (1) 

Such an equation is without any doubt among the most investigated ones in 
the qualitative theory of ordinary differential equations, for its applications in 
mechanics and electric circuits theory. It may be considered the starting point 
in the theory of limit cycles, becoming sort of a benchmark for new methods. 
Its relationships with different classes of systems, as quadratic ones (which can 
be reduced to suitable Lienard systems) made it even more relevant in the study 
of planar systems. 

We do not intend to make a survey on this widely investigated area. Our 
attempt is rather directed to focus on the results of Giovanni Sansone and to 
show how in some sense those results were giving a crucial contribute for the 
development of this field. After such equation was proposed by Lienard in 1928 
[TU], who proved the first uniqueness theorem for limit cycles, it was necessary 
to wait until 1942 to find another significant result in this direction, due to 
Levinson and Smith [9]. But it was only with the contribution of Sansone, 
starting from his celebrated paper of 1949 [TB], that the relevance of this problem 
was recognized, not only from the Italian school. We just recall the famous 
paper of Filippov (6j, which was devoted to the problem of existence of limit 
cycles and therefore will be not examined in this context. For the enormous 
contribution of the Chinese School, which for long time was not known out 
of China for language problems, we refer to the monograph of Ye Yan Qian 
"Theory of limit cycles" [H]. Finally, a survey updated until the 70's is due to 
Staude PO*. Coming back to the work of Sansone, we believe that our, in some 



sense historical, approach will be of some interest, because it shows how deep 
results may come out from simple geometrical ideas in a far to be easy field of 
research as is the problem of uniqueness for periodic solutions. In particular, 
we would like to describe the three approaches applied by Sansone in such a 
field: the energy integral, the divergence integral and, implicitly, the rotation of 
the vector field along rays. The first one was not new, since it appeared in the 
very first paper about Lienard equation [TU] , which gave to such an equation its 
name. The second one was probably used for the first time just by Sansone, as 
well as the third one, even if this approach is usually credited to Massera for 
his geometrical idea centered on the rotation of the vector field on rays. 

In section 2 we give an account of the evolution of the theory starting from 
Levinson-Smith paper. Some of the more recent developments in this direction 
will be also presented. We know that not all the authors which deserve to be 
quoted will appear and we apologize for this. 

In section 3 we discuss some problems arising from the geometrical approach 
to uniqueness. We also present a new result in the light of the ideas of Sansone 
and Massera. Such a result can be applied to a class of planar differential 
systems, not equivalent to second order differential equations. This extends a 
recent result by Ciambellotti also based on Massera's geometrical approach [4]. 

In section 4 we point out that the field still uncovered by existing uniqueness 
theorems contains strikingly simple, but non-trivial equations. Lins-De Melo- 
Pugh's conjecture, still open for n > 3, resists to researchers' efforts since 1976. 



2 Uniqueness theorems 

Let us consider a Lienard equation in the phase plane, 

i = y, y = -9{x) - f{x)y. (2) 
We shall also consider the usual system, known as Lienard plane, 

x^y-F{x), y = -g{x), (3) 

where F{x) = f{s)ds. Throughout all of this paper we assume f{x) to 
be continuous and g{x) to be locally lipschitzian on their common domain, so 
to have existence and uniqueness of solutions. In the following we also write 
G{x) = /q g{s)ds, referring to 

E{x,y) = Gix) + ^ 

as the energy function of both systems ([2| and (|3| . Some of the basic properties 
and techniques concerned with limit cycles' existence and uniqueness appear 
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already in Lienard's paper [TU], which may be considered a milestone in this 
field, but after this paper, the first relevant result concerning the uniqueness 
of limit cycles is due to Levinson and Smith [9]. In particular, the authors 
presented the following results. We write F{+oo) for °° f{x)dx. Similarly 
for G. 

Theorem 1. Assume that xg{x) > for x ^ 0. If 

• there exist 5^ < < such that G{S^) = G{5+); 

• f{x) < m {5-,5+); 

• fix) >0 in (-00, U [6+ , +oo); 

• G(±oo) = F(+oo) = +oo; 

then the system ^ has exactly one limit cycle. 

Such a theorem is obtained as a special case of a more general one concerning 
the system 

x = y, y ^ -g{x) ~ f{x,y)y. (4) 

Restricting to classical Lienard systems, the same paper contains also a result 
strongly depending on symmetry properties of f{x) and g{x). Even if it is not 
explicitly stated as a corollary, we can summarize the containt of section 4 in 
[S] in the following statement. 

Theorem 2. Assume f{x) even, g{x) odd, and xg[x) > for x ^ 0. If 

• there exists xq > such that F{x) < in (0, xq); 

• F{x) > and increasing in (xo,+oo); 

• G(+oo) = F(+oo) = +oo; 

then ^ has exactly one limit cycle. 

The proof did not introduce relevant novelties, being based on Lienard's 
approach. In fact, theorem ([2|, under the condition g{x) = x, was proved 
by Lienard in his celebrated paper [TD]. Both Lienard's and Levinson-Smith's 

f)ofs are based on the fact that the energy's derivative along the solutions of 
is E{x,y) = —F{x)g{x), so that the integral of —F{x)g(x) along a cycle is 
zero. Splitting the integral in several steps, where the integration is performed 
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with respect to x or y, and comparing analogous sub-integrals of distinct cycles, 
one proves that if 71 and 72 are two concentric cycles, with periods Ti and 



so obtaining a contradiction. 

A few years after Levinson and Smith, Sansone proved the following theorem, 
where the symmetry assumption on F{x) was replaced by a weaker one |16) . 

Theorem 3. If g{x) = x, and 

• there exist 5^ < Q < 5^ such that f{x) < for x G (S^ , S^); 

• either f{x) > in (5+, +00) or f{x) > in (—00, 6^); 

• there exists A > such that F{A) = F(-A) = 0; 

• F(+oo) = +00, or F(— 00) = —00; 

then the system ^ has exactly one limit cycle, which is stable. 

In the same direction, for more general classes of systems, we can find the 
recent results of Xiao and Zhang Zhi-fen [22], [23], Carletti and Villari [2], 
Sabatini and Villari [TS], Carletti [11. 

Next theorem was based on a different principle |16j . 
Theorem 4. If g{x) = x, and 

• there exists 6 > such that f{x) < for x G {—S,S), f(x) > in 
{-00, -S) U {6, +00); 

• -F(+oo) = +00, or F(—oo) = —00; 

then the system ^ has exactly one limit cycle, which is stable. 

The above theorem, even if a special case of theorem [T] introduced a sta- 
bility argument often re-used by several authors in successive papers. Sansone 
proves that if a T-periodic limit cycle "f{t) exists, then it is attractive, since the 
divergence integral div {'j{t))dt is negative. Then, since two adjacent con- 
centric limit cycles cannot be both attractive, the uniqueness follows. Such a 
theorem, even if not explicitly stated, proves also the cycle's hyperbolicity, that 
is an important feature in relation to perturbation problems. 
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An equally innovative result obtained by Sansone on limit cycle's uniqueness 
is virtually unknown to most researchers. It was exposed in a talk, together some 
other results about Lienard equation, whose containt appeared in [17], section 
2. We report here the uniqueness theorem for Lienard equation. 

Theorem 5. // g{x) — x, and 

• there exist < < such that f{x) < for x G {5^ , 

• either f{x) > in (5+, +oo) or f{x) > in (— oo, S^); 

• fiS-)^fiS+)^0; 

• 1/(^)1 < 2; 

• f{x) is non-increasing in {—oo,6^) and non- decreasing in ((5+,+oo) ; 

then the system has exactly one limit cycle. 

The proof was based first on the trasformation of the system ([3| into polar 
coordinates (r, 6), then on the study of the function 

^ ^ ^(rcosejjAr^ 

dd ^ ' 1 + f{r cose) sm9cos9 

A generalization of Sansone's result is due to Conti, who used the map 

{x,y) {sign{x)^/G{x),y) 

known as Conti-Filippov transformation, to transform the system ([2]) or ([s]) into 
a new system with g{x) replaced by x. 

The weakness of last Sansone's result is the assumption |/(a;)| < 2 (which 
avoids the denominator's vanishing), clearly not satisfied by polynomials. In 
Sansone's approach it was not clear whether such an obstacle could be overcome. 
Such an assumption was removed by Massera, when he was visiting Florence 
and working with Sansone. Actually, Massera's main contribution consisted in 
turning Sansone's analytical approach into a geometrical one. Massera observed 
that the monotonicity conditions on f{x) imply that the vector field rotates 
clockwise, as r increases, along rays contained in the half-plane a; > 0, counter- 
clockwise, as r increases, along rays contained in the half-plane a; < 0. If 7 is 
star-shaped, then ever curve 7" obtained from 7 by a K-dilation (homothety) is 
as well star-shaped. The family 7^ is a foliation of the punctured plane, and can 
be assimilated to the family of level curves of a Liapunov function for 7. The 
vector field rotation forces the orbits of ([2| to enter such curves, approaching 7. 
This in turn allows to show that every limit cycle is attractive, so preventing the 
co-existence of concentric limit cycles, by Sansone's attractivity argument. On 
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the other hand, if the hmit cycle is not star-shaped, the curves 7fc can intersect 
with each other, and for values of k close to 1 they even intersect the limit cycle 
7, so that proving 7's attractivity is not immediate. It is amazing to observe 
that in the original paper Massera deals with this point by just writing II faut 
remarquer que le raisonnement resterai valable meme si 7 n'est pas etoile par 
rapport a I'origine: dans ce cas 7 et 7^ pourraient avoir des intersections. 

It is worth noticing that Sansone's condition |/(a;)| < 2 implies that for every 
non-trivial orbit the angular speed does not vanish. In particular, this occurs 
at cycles, which hence have to be star-shaped. 

It is also surprising to observe that the cycle's star-shaped-ness could have 
been easily proved by using the monotonicity conditions on f{x), as shown in 
pS] . p. 225. A more geometrical proof of the cycle's star-shaped-ness was 
developed by Villari, and appeared for the first time in [3]. A second objection 
to Massera's geometrical approach is that, in general, the vector field could 
rotate up to reach a position opposite to the original one (tangent to 7) so that 
the successive position would take the orbits to leave the external 7fc's. This 
cannot occur for second order equations, as will be shown in section 3. 

We emphasize that the cycle's star-shaped-ness does not imply that the 
angular speed be of constant sign on every orbit. In fact, even for the simplest 
example of Lienard equation with a limit cycle, i.e. Van der Pol's equation. 



infinitely many orbits passing through the second and fourth orthants change 
angular speed. We may consider the information on the cycle's star-shaped- 
ness also as information about the cycle's location: it is contained in the region 
x'^ — xyf{x)+y'^ > 0. Other extensions were given by Carletti, Rosati and Villari 
[3] and Ciambellotti [1] , but in both papers the geometrical idea of Massera is 
not really modified. 

At the end of this review, we cite a result that renews the divergence ap- 
proach in a non-traditional way. In some recent papers new classes of so-called 
stability operators were introduced, in order to study cycles' hyperbolicity. Such 
operators are just functions "^{x,y) with the property that if 7 is a T-periodic 
cycle and 



then 7 is hyperbolic, the stability character being given by the sign of the above 
integral. The divergence is clearly one of such operators. Other ones are the 
curvature of the orthogonal vector field, or the function 



X + e(x^ — l)x -t- a; = 0, 




[V, W]AV 
V AW 



V AW ^0, 
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where [V, W] is the Lie bracket of V and W, and V AW is the determinant of 
the matrix having V and W as rows ([7], [5]). The function v has the foUowing 
relevant property 



(div V){j{t)) dt = ^ ' (7(0) dt, 



T 

In particular, if < in a domain (not identically vanishing on any cycle), then 
such a domain cannot contain two concentric limit cycles, since they both should 
be attractive. This is just the original argument by Sansone, applied replacing 
the divergence with v. We emphasize that using v allows also to prove the cycle's 
hyperbolicity, which is not a consequence of Sansone-Massera's approach. 

Finally, we observe that Massera theorem is a special case of such a result, 
obtained by taking W{x,y) = {x,y). In fact, as shown in Corollary 6 of [5], 
given an arbitrary planar differential system 

x^P{x,y), y^Q{x,y), (5) 

taking W{x, y) = {x, y), gives, for Q ^0, 

^^ P{xQ,+yQy)-QixP,+yPy) ^ "-{^p) ^ r / d y\ 
yP — xQ yP — xQ yx — xy \dr x ) 

In other words, if v does not change sign, then the vector field rotates mono- 
tonically along rays as r increases. Moreover, if 

X = P{x,y) =y, y^Q{x,y) = -x-yf{x), 

then ^ 

^ x/''(x) ^ (7) 

+ xyf{x) + 

Hence, if the denominator x'^ + xyf{x) + y^ does not vanish, j/'s sign is the same 
as that of xf'(x), rc- finding Sansone-Massera's condition. 



3 A new result 



Let us call radial angular monotonicity, RAM, the property of the vector field 
to rotate monotonically along rays as r increases. An implicit assumption of 
Massera-like theorems is the fact that limit cycles rotate clockwise around the 
origin. This is a consequence of i — y, which forces cycles cross the y-axis 
rotating clockwise. If RAM holds, this gives the attractivity of every cycle, 
hence its uniqueness. In fact, the same assumption on a cycle rotating counter- 
clockwise generates repulsiveness, rather than attractiveness. This allows as 
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well to prove uniqueness, if RAM holds. On the other hand, the possible co- 
existence of limit cycles rotating both clockwise and counter-clockwise does not 
allow to prove uniqueness under the only RAM hypothesis. Counterexamples 
with finitely many and infinitely many cycles have been given in |14) . We 
report here a polynomial system satisfying RAM, with two limit cycles rotating 
in different ways. 




Such a system has two star-shaped limit cycles coinciding with the circles 
y2 _ 3->/5 ^2 y2 _ 34V5 ^ rpj^g internal one is an attractor, the external 
one is a repellor. The vector field rotates clockwise along every ray (see figure 
!)• 



Figure 1: The system ([s]) has two limit cycles. 
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Hence, in order to prove a Massera-like theorem valid for systems, one has 
to assume some additional hypothesis that prevents the coexistence of opposite 
rotations. This docs not necessarily have to hold in all of the plane. It is 
sufficient, to prove that it occurs in a region containing all cycles, or, as in 
[131 m], in the region bounded by two adjacent limit cycles. In [T31 [T3], this 
was proved for systems of the type 

i = 2/, y ^ -x-y(j){x,y) 

under the assumtpion that 0's level sets be star-shaped. This does not imply 
the orbits' star-shaped-ness. 

We propose here some simple additional conditions that allow to prove 
uniqueness if RAM holds. We state next theorem for a system defined on 
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the whole plane, but it can be easily adapted to an arbitrary star-shaped set 
containing the origin. We assume existence and uniqueness of solutions to ([s]). 
Let us set 

a{x, y) = P{xQ.^ + yQy) - Q{xP^ + yPy). (9) 
In next theorem, the interval (a, h) may be a half-line or even all of IR. 

Theorem 6. Assume that ^ to have a unique equilibrium point at the origin. 
Assume a{x, y) > in fl. If one of the following conditions holds 

(1) for all {x, y) G IR^, {x, y) ^ (0, 0), for all r > 0, the vectors V{tx, ry) and 
V{x,y) are linearly independent; 

(2) there exists a , open Jordan curve curve rj : (a, 6) — >■ IR^ such that 
lims^<j+ |77(s)| = 0, lim^^f,- |77(s)| = +oo, and 9{t]{s)) > (< 0); 

then the system ^ has at most one limit cycle in fl. 

Proof. According to the sign of a is the sign of the radial derivative of 

y 

— . Working as in :3| , one proves that every cycle is star-shaped The existence of 

a single equilibrium point implies that all cycles are concentric. Then we work 
as follows. 

Point (1). It is a modification of the original argument by Massera. If a 
cycle 7 exists, as in [12j we may consider the family of curves 7^ obtained from 
7 by means of K-dilations (homotheties). In particular, 7 = 71. Without loss of 
generality, we may assum 71 to rotate clockwise. Then, the vector field's rotation 
along rays implies that, for small positive values of k — 1, V points towards the 
interior of 7^. Assume, by absurd, the existence of k* and {x*,y*) such that 
V{x* , y*) points outwards 7^;* . Then, by continuity, there exists and (a;^, y^) 
such that V{x^ ,y^) is tangent to 7„+. The rotation's monotonicity along rays 
implies that this occurs first at a point y+) where V{x^ , t/+) is parallel and 
opposite to V{xi,yi), with {xi,yi) G 7. This contradicts the hypothesis (1). 
Internal attractivity can be proved in the same way. 

Point (2). Assume, by absurd, the existence of two distinct cycles 7i, 7e, 
with 7j encircled by je- Without loss of generality, we may assume them to be 
adjacent. Since both are star-shaped, 9 does not change sign on any of them. 
Moreover, since they both cross the curve rj, both have to contain a point where 
^ > 0. Hence e{-fi) > and e{je) > 0. The condition RAM implies that both 
are attractive, which is a contradiction. Jjt 

Both hypotheses (1) and (2) are satisfied by Lienard systems. Point (1) since 
i = y on the curve 77(5) = (0, s), defined on (0, +00). 
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Next corollary extend Ciambellotti's result Consider the following family 
of systems 



X = Hv), y = -f{x)l{y) - 2^ hj{x)mj{y). 

3=0 



(10) 



Corollary 1. Assume k{y), l{y), to he d-homogeneous functions, hj{x) j- 
homogeneous, mj{y) {d — j) -homogeneous, for j — 1, . . . ,7i, f{x) differentiable. 
If yk{y) > and yl{y) > for y ^ 0, xf'{x) > for x 0, then the system 
[Tm has at most one limit cycle. 



Proof. The system ( 10 ) satisfies the hypothesis (2) of theorem |6j taking 
77(5) = (0, s), s € (0, +00). Let us compute the function a for the system ([To|: 



a = k 



fl + J2hj 



3=0 



yk' 



cf'lk — k xh'jTUj — fyl'k — k hjym'^ + flyk' + yk' hj 



3=0 



= —xf'lk — k jhjnij — dflk — k hj {d — j)mj + dflk + dk h 



3=0 



-xf'lk + k 



3=0 



-xf'lk. 



The sign of —xf'{x)l{y)k{y) is the same as that of -~xf'(x), since both k{y) and 
l{y) have the sign of y. Hence the sign of a is that of —xf'{x), as in formula 
0. * 



An example is given by the the system 

i = y^ y = (5a;2 - l)z/ - a;^ - xy^ , 
illustrated in figure 2. 



(11) 
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Figure 2: The system (11) has a unique hmit cycle 




4 An open problem 



At the end of this paper, after having hsted several theorems of seemingly great 
generality, we would like to point out a very simple problem still waiting for a 
solution. Let us consider the Lienard equation with /(x) cubic polynomial, 



X + [ax^ + bx^ + cx + d)x + a; = 0, 



a, b,c,d <E ]R. 



Such an equation is clearly out of the reach of any of the above results. Proving 
the limit cycle's uniqueness for such a class of equations would give an answer 
to Lins-De Melo-Pugh's Conjecture for n = 3 [11]. As it is well known, such a 
conjecture is part of the famous XVI Hilbert's problem, re-proposed by Smale 
as one of the Mathematical problems for the next century |19| . The conjecture 
for n = 2 was proved by Lins-De Melo-Pugh in their paper, but they were not 
aware of the fact that this problem was actually already solved by Zhang-Zhifen 
in 1958. This is a very deep result, but the for long time the proof was only in 
Chinese, until an English version appeared [23]. 
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